Novel Approach to Real Polynomial Root-finding and Matrix Eigen-solving
نویسنده
چکیده
Univariate polynomial root-finding is both classical and important for modern computing. Frequently one seeks just the real roots of a polynomial with real coefficients. They can be approximated at a low computational cost if the polynomial has no nonreal roots, but typically nonreal roots are much more numerous than the real ones. We dramatically accelerate the known algorithms in this case by exploiting the correlation between the computations with matrices and polynomials, extending the techniques of the matrix sign iteration, and exploiting the structure of the companion matrix of the input polynomial. We extend some of the proposed techniques to the approximation of the real eigenvalues of a real nonsymmetric matrix.
منابع مشابه
New progress in real and complex polynomial root-finding
Matrix methods are increasingly popular for polynomial root-finding. The idea is to approximate the roots as the eigenvalues of the companion or generalized companion matrix associated with an input polynomial. The algorithms also solve secular equation. QR algorithm is the most customary method for eigen-solving, but we explore the inverse Rayleigh quotient iteration instead, which turns out t...
متن کاملTR-2006007: Root-Finding with Eigen-Solving
We survey and extend the recent progress in polynomial root-finding based on reducing the problem to eigen-solving for highly structured generalized companion matrices. In particular we cover the selection of the eigen-solvers and the matrices and the resulting benefits based on exploiting matrix structure. No good estimates for global convergence of the basic eigen-solvers have been proved, bu...
متن کاملRandomized Preprocessing of Homogeneous Linear Systems of Equations
Our randomized preprocessing enables pivoting-free and orthogonalization-free solution of homogeneous linear systems of equations, which leads to significant acceleration of the known algorithms in the cases of both general and structured input matrices. E.g., in the case of Toeplitz inputs, we decrease the estimated solution time from quadratic to nearly linear, and our tests show dramatic dec...
متن کاملTR-2009008: Randomized Preprocessing of Homogeneous Linear Systems of Equations
Our randomized preprocessing enables pivoting-free and orthogonalization-free solution of homogeneous linear systems of equations, which leads to significant acceleration of the known algorithms in the cases of both general and structured input matrices. E.g., in the case of Toeplitz inputs, we decrease the estimated solution time from quadratic to nearly linear, and our tests show dramatic dec...
متن کاملTR-2009014: Randomized Preprocessing of Homogeneous Linear Systems
Our randomized preprocessing enables pivoting-free and orthogonalization-free solution of homogeneous linear systems of equations. In the case of Toeplitz inputs, we decrease the solution time from quadratic to nearly linear, and our tests show dramatic decrease of the CPU time as well. We prove numerical stability of our randomized algorithms and extend our approach to solving nonsingular line...
متن کامل